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ABSTRACT 

The recent discoveries of massive planets on ultra-wide orbits orbiting HR 8799 (Marois et al. 
2008) and Fomalhaut (Kalas et al. 2008) present a new challenge for planet formation theorists. 
Our goal is to figure out which of three giant planet formation mechanisms — core accretion (with or 
without migration), scattering from the inner disk, or gravitational instability — could be responsible 
for Fomalhaut b, HR 8799 b, c and d, and similar planets discovered in the future. This paper presents 
the results of numerical experiments comparing the long-period planet formation efficiency of each 
possible mechanism in model A star, G star and M star disks. 

First, a simple core accretion simulation shows that planet cores forming beyond 35 AU cannot 
reach critical mass, even under the most favorable conditions one can construct. Second, a set of 
N-body simulations demonstrates that planet-planet scattering does not create stable, wide-orbit 
systems such as HR 8799. Finally, a linear stability analysis verifies previous work showing that 
global spiral instabilities naturally arise in high-mass disks. We conclude that massive gas giants 
on stable orbits with semimajor axes a > 35 AU form by gravitational instability in the disk. We 
recommend that observers examine the planet detection rate as a function of stellar age, controlling 
for the planets' dimming with time. Any age trend would indicate that planets on wide orbits are 
transient relics of scattering from the inner disk. If planet detection rate is found to be independent 
of stellar age, it would confirm our prediction that gravitational instability is the dominant mode of 
producing detectable planets on wide orbits. We also predict that the occurrence ratio of long-period 
to short-period gas giants should be highest for M dwarfs due to the inefficiency of core accretion and 
the expected small fragment mass (~ 10Mj up ) in their disks. 

Subject headings: planetary systems — planetary systems: formation — accretion disks — stars: 
formation — instabilities 
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1. INTRODUCTION 

The first direct images of extrasolar gas giants present 
planetary scientists with a new mandate: construct mas- 
sive planets on extremely wide orbits. Before November 
2008, when the existence of planets surrounding Fomal- 
haut (Kalas et al. 2008) and HR 8799 (Marois et al. 2008) 
was revealed, theorists had reached a near-consensus that 
the core accretion mechanism could explain nearly all ob- 
served properties of planets discovered by radial veloci- 
ties. Simple core accretion-based models reproduce the 
planet-metallicity correlation (Fischer & Valenti 2005) 
and the planet-silicon correlation (Robinson et al. 2006), 
and the rarity of giant planets orbiting M dwarfs (Laugh- 
lin et al. 2004; Johnson et al. 2007). Unfortunately, with 
the exception of Ida & Lin (2004), these models treat 
only planets forming near the ice line and ignore planet 
formation beyond ~ 5 AU — for the simple reason that 
the vast majority of known radial-velocity planets orbit 
within 3 AU of their host stars. Estimates of the planet 
occurrence rate beyond 5 AU (e.g. Cumming et al. 2008) 
rely on extrapolation of the statistics on shorter-period 
orbits. 

One can immediately foresee problems applying the 
core accretion model — in which a massive solid core 
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destabilizes the surrounding disk gas to accrete an 
atmosphere — to Fomalhaut b, located 119 AU from the 
host star, and HR 8799 b, c and d, with projected stellar 
separations of 24, 38 and 68 AU. All of these new plan- 
ets are unequivocally massive and gaseous: the smallest 
is Fomalhaut b with M < 3Mj up . Even forming the 
small and solid-rich Uranus and Neptune (14.2 and 17.5 
M ffl , respectively) in situ is difficult: Levison & Stew- 
art (2001) show that the planets cannot form in their 
current orbits at 19 and 30 AU because the accretion 
timescale is far longer than the planetesimal excitation 
timescale. If accretion timescale prohibits medium-mass 
planet formation at 19 AU, the formation of a gas giant 
at > 100 AU by core accretion hardly seems possible. 

There are a few ways to increase the efficiency of core 
accretion. Perhaps the protoplanetary disks surround- 
ing HR 8799 and Fomalhaut held far higher planetesimal 
masses than the solar nebula, and were thus better able 
to form solid cores at large radii. Since both stars host 
massive, bright debris belts (Kalas et al. 2005; Chen et 
al. 2009), this hypothesis requires investigation. It's also 
possible that these relatively massive stars (~ 1.5M Q , 
~ 2.0M Q respectively) decreased the protoplanet orbital 
timescale in the outer disk enough for core accretion to 
proceed. Perhaps it is not a coincidence that the only 
directly imaged extrasolar planets orbit A stars. In this 
work we use model A, G and M-star disks to test the 
possibility that core accretion at large radii could work 
for HR 8799 and Fomalhaut where it fails for lower-mass 
stars. 

What are the other possible formation mechanisms for 
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gas giants on wide orbits? For truly massive planets 
and in the case of HR 8799, the best-fit masses of 7 Mj up , 
10 Mj up and 10 Mj up are near the Y dwarf minimum 
mass (Zuckerman & Song 2009) — the simplest possibility 
is fragmentation of the collapsing protostellar core, anal- 
ogous to binary star formation (e.g. Cha & Whitworth 
2003). However, several lines of observational evidence 
show that HR 8799 and Fomalhaut are not simply binary 
or multiple stars with high mass ratios: 

1. Resolved images show that both stars host debris 
disks that are confined by the planet (s), indicating 
that the planets and debris disks are coplanar relics 
of protostellar disks (Kalas et al. 2005; Chen et al. 
2009; Su et al. 2009). 

2. Spectroastrometric measurements of the spin-orbit 
alignment of Fomalhaut show that the star's rota- 
tion axis is perpendicular to the debris disk and 
planet orbit, further indicating that Fomalhaut b 
formed in a disk (Le Bouquin et al. 2009). The ro- 
tational velocity of HR 8799 indicates that it, too, 
is nearly pole-on and perpendicular to the planet 
orbits (Reidemeister et al. 2009). 

3. The three companions to HR 8799 orbit the central 
star, rather than forming a double-double system 
consistent with cloud core fragmentation (Clark & 
Bonnell 2005). 

As molecular cloud core fragmentation is unlikely, we 
will treat the companions to HR 8799 and Fomalhaut as 
bona fide planets (whether they burn deuterium or not) 
and search for formation mechanisms that take place in 
protoplanetary disks. The two remaining possibilities 
are planet-planet scattering, in which planets form in 
the inner disk and are ejected to wide orbits by a yet- 
unseen inner planet, and gravitational instability, where 
the protostellar disk fragments into bound clumps that 
eventually contract to near Jupiter's radius. Here, too, 
we run into problems. Planet scattering tends to cre- 
ate either unstable systems or extremely eccentric orbits 
(Raymond et al. 2008; Veras et al. 2009), whereas the 
orbits of HR 8799 b, c and d and Fomalhaut b appear to 
be nearly circular. There is considerable debate in the 
literature over whether disk fragments formed in gravita- 
tional instabilities survive to become protoplanets or are 
destroyed by accretion onto the star or background shear 
flows (Durisen et al. 2007; Vorobyov & Basu 2006); fur- 
thermore, a previous investigation of planet formation 
on wide orbits by Boss (2006) noted that spiral insta- 
bilities may wind too tightly at large radii for the disk 
to fragment at all. Nonetheless, two points in favor of 
gravitational instability are that spiral waves — whether 
they fragment the disk or not — develop most easily in 
the outer parts disks, and that when fragments do form 
they are often near 10Mj up (Boley 2009; Stamatellos & 
Whitworth 2009). Gravitational instability is therefore 
perfectly positioned to create massive planets with ex- 
tremely long-period orbits. 

In this paper we present numerical experiments that 
test the viability of core accretion, planet-planet scatter- 
ing and gravitational instability. The paper is organized 
as follows. In ^2) we describe the disk models that form 
the foundation of our numerical experiments. In $3l we 



use a simple core-accretion simulation to show that mas- 
sive core growth on ultra- wide orbits is not possible. We 
present N-body simulations assessing the likelihood of 
scattered planets landing in stable, wide orbits in $4] The 
disk stability analysis in <J5] confirms that our maximum- 
mass model disks can become unstable to global two- 
armed spiral modes. Finally, we present our conclusions, 
along with a proposed observational experiment to con- 
firm that gravitational instability is the dominant mode 
of planet formation on wide orbits, in Sj6] 

2. DISK MODEL 

We begin by describing the simple disk model that 
forms the basis for our three numerical experiments. In 
order to facilitate comparisons of planet formation effi- 
ciency among different types of stars, we have deliber- 
ately chosen a simple, passively illuminated disk whose 
temperature profile depends only on the temperature and 
luminosity of the host protostar. 

We use the method of Chiang & Goldreich (1997) to 
construct flared protostellar disks that are in radiative 
equilibrium with their host protostars. In the Chiang & 
Goldreich models, dust in a superheated surface layer 
reprocesses the incident stellar radiation, re-radiating 
~ 1/2 of the absorbed flux onto the disk interior at in- 
frared wavelengths. The midplane temperature of a disk 
optically thick to visible stellar radiation is 
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In Equation [TJ a is the distance from the star, is the 
protostar radius, T, is the effective temperature and a is 
the disk flaring angle: 
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The first term in Equation[2]can be neglected in the limit 
a > R*. 

The visible photosphere height is given by 
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where h is the pressure scale height. The quantity T c 
defines the escape temperature of the disk gas from the 
stellar surface: 

In EquationJU Af* is the star mass, /j, is the mean molec- 
ular weight of the nebula gas in grams and k is Boltz- 
mann's constant. We take H/h — 4 as a constant among 
all model disks. Best-fit SEDs of protostellar disks have 
values H/h between 2 and 5 (Chiang et al. 2001). 

We test the planet formation efficiency in disks sur- 
rounding the types of stars known to host planets (spec- 
tral types A-M) by defining three representative proto- 
stars: an A star of mass 1.5M Q , a G star of mass I.OMq, 
and an M star of mass 0.5M Q . The remaining inputs to 
Equations[THH the star radius and effective temperature, 
are taken from published pre-main-sequence evolution- 
ary tracks (D'Antona & Mazzitelli 1994). We choose a 
stellar age of 10 5 years as the onset of planet formation. 



Planets on wide orbits 



3 



Our stars are older than the 100 km planetesimal forma- 
tion timescale (Barnes et al. 2009) yet still young enough 
to harbor massive disks, which improves the chances of 
planet formation. Figure [1] (top panel) shows the mid- 
plane temperature as a function of radius for each pro- 
tostar mass. 

We must now calculate the mass available for planet 
formation. The Chiang & Goldreich temperature profile 
is independent of the disk surface density as long as the 
disk is optically thick. We therefore use the constraint 
of local gravitational stability to set the surface density 
profile. The Toomre Q parameter, 



(5) 



ttGE' 



measures the stability of the disk to axisymmetric per- 
turbations. In Equation [5l k is the epicyclic frequency, 
c s is the midplane sound speed and X is the gas surface 
density. Rings of enhanced surface density, with cor- 
responding rarefactions, will grow exponentially in any 
part of the disk with Q < 1. With c s determined by 
the radiative equilibrium model and the assumption of a 
nearly Keplerian disk, we can recover the surface density 
profile from Equations Q] [5] by specifying a value of Q: 
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Since planets on wide orbits are known to be difficult 
to build (Kenyon & Bromley 2008; Whitworth & Sta- 
matellos 2006), we construct our disks with the follow- 
ing question in mind: Are there any circumstances, how- 
ever optimistic, under which massive planets between 20 
and 120 AU can form and survive? We therefore outfit 
each star with a "maximum-mass nebula" , a disk that 
contains as much mass as it can plausibly hold while 
maintaining the disklike structure that is critical for core 
accretion. The maximum-mass nebulae have Q = 1.5 
everywhere in the disk. 

Furthermore, since most known planets are not on wide 
orbits and do not require extremely massive disks to 
form, we also construct "medium-mass nebulae" with 
Q = 8 for each star. The medium- mass nebulae allow 
us to verify that our calculations of planet formation ef- 
ficiency are consistent with previous studies of both core 
accretion and gravitational instability (e.g. Laughlin et 
al. 2004; Kenyon & Bromley 2008; Kratter et al. 2008; 
Stamatellos & Whitworth 2009). 

The final step in constructing our disk models is 
to account for the viscous accretion that happens be- 
fore the onset of planet formation, which tends to flat- 
ten a strongly centrally peaked surface density profile 
(Lynden-Bell & Pringle 1974; Dodson-Robinson et al. 
2009). We calculate the viscous diffusion timescale, over 
which surface density decreases by a factor of ~ 2, as 
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where a v isc is the viscous efficiency in the standard a-disk 
model (Shakura & Syunyaev 1973). Following Dodson- 
Robinson et al. (2009) and Lyra et al. (2008), we take 
«visc = 0.002. Our final surface density is determined by 



reducing the fiducial value according to 
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where idiff is given in years. Our method of computing 
the diffusion timescale includes the assumption that vis- 
cous torques provide some angular momentum transport, 
but stellar radiation dominates the disk energy budget. 
The effect of diffusion is most profound in the inner disk 
where idiff is smallest. The slightly reduced surface den- 
sity in the inner disk provides some additional gravita- 
tional stability to the disk as a whole. 

FigureQ] (bottom panel) shows the solid surface density 
profiles of the maximum-mass and medium-mass nebu- 
lae for each star. Each disk has an outer boundary of 
100 AU. Gas/solid ratios are taken from the disk molec- 
ular inventory at solar composition of Dodson-Robinson 
et al. (2009). The maximum-mass disks hold between 1/4 
and 1/3 of the system mass and may undergo global grav- 
itational instabilities that aid planet/planetesimal forma- 
tion (e.g. Adams et al. 1989; Rice et al. 2004; Lodato & 
Rice 2005), which we examine in <J5] Such massive disks 
would be either in the middle or just at the end of their 
epoch of accretion from the protostellar clump, consis- 
tent with our assumed stellar age of 10 5 years. Although 
we do not examine accretion-triggered instability, it is 
important to note that the accretion rate onto the disk 
can be a determining factor for fragmentation at large 
radii (Vorobyov & Basu 2006; Boley 2009; Clarke 2009; 
Rafikov 2009). The medium-mass nebulae hold < 10% 
of the total system mass, consistent with disk masses de- 
rived from sub-mm observations (Andrews & Williams 
2007). 

Having created our disk models, we now test the suit- 
ability of each for forming planets. 

3. EXPERIMENT 1: CORE ACCRETION 

Since the astronomical community has reached a near- 
consensus that core accretion is the dominant mode of 
planet formation — explaining such observables as the 
planet-metallicity correlation (Gonzalez 1998; Fischer & 
Valenti 2005), the planet-silicon correlation (Robinson et 
al. 2006), and the paucity of gas giants orbiting M dwarfs 
(e.g. Johnson et al. 2007), we begin our investigation of 
planets on wide orbits with an assessment of the ability 
to form critical-mass cores at large radii. The quantity of 
interest is the maximum possible protoplanet core mass 
as a function of distance from the star. If the core can 
reach 10M®, the canonical critical value for destabilizing 
the surrounding nebula to accrete a massive atmosphere 
(Mizuno 1980), while the gas disk is still present, it has a 
chance of forming a gas giant. We adopt a disk lifetime of 
5 Myr, the maximum star cluster age with an apprecia- 
ble disk fraction according to Currie et al. (2009). Here 
we are being generous: most disks have far shorter ob- 
served lifetimes, but we must account for the possibility 
that planets on wide orbits are relics of long-lived disks. 

We use the feeding zone approximation (Safronov 1969; 
Lissauer 1993), in which a solitary core accretes planetes- 
imals at the rate 
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where M is the core mass, R is the core radius, £ is the 
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surface density of planetesimals, is the Keplerian fre- 
quency and (v) is the RMS planetesimal velocity. Assum- 
ing solitary cores tends to increase the maximum planet 
core mass, since oligarchic growth rates are much slower 
than those of a single embryo in a swarm of planetes- 
imals. The feeding zone extends 4 Hill radii on cither 
side of the protoplanet (Kary & Lissauer 1994), where 
the Hill radius Rh demarcates the protoplanet 's Roche 
lobe in the spherical approximation: 

/ M \ 1/3 
RH - a {3Mj ■ (10) 

In Equation [lOl a is the planet's semimajor axis and M* 
is the star mass. We begin the calculation with a 0.1M® 
seed core and use a fourth-order Runge-Kutta method to 
integrate Equation [9] for 5 Myr, our adopted maximum 
disk lifetime. We consider only the solid core growth and 
do not treat gas accretion. 

To calculate the gravitational focusing factor (the ratio 
of the core escape velocity to the mean planetesimal ve- 
locity), we apply the simple assumption that (v) = RrQ, 
so that planetesimals are at the Hill velocity. Accre- 
tion is therefore on the threshold between dispersion- 
dominated, where only two-body forces (protoplanet 
and planetesimal) are important, and shear-dominated, 
where the entire star-planet-planetesimal system must be 
treated (Greenzweig & Lissauer 1992). Since three-body 
accretion rates tend to be lower than two-body rates, our 
core accretion model may provide optimistically high es- 
timates of the core accretion efficiency. However, as we 
shall see, core accretion cannot form gas giants at large 
radii even under the extremely favorable conditions we 
have constructed here. 

Figure [2] shows the maximum planet core mass as a 
function of distance from the star. In the medium-mass 
A and G star nebulae, giant planet formation is com- 
fortably confined to a narrow belt near 10 AU: there 
is no chance of forming a super-Jovian planet beyond 
20 AU. The maximum core mass barely reaches the crit- 
ical lOAf^ in the G star disk, confirming that giant planet 
formation by core accretion is a threshold phenomenon. 
A slight downward adjustment in disk mass, metallicity 
or lifetime would be enough to prevent any giant planets 
from forming — and indeed, the estimated planet occur- 
rence rate between and 5 AU is only 10% for solar-mass 
stars (Cumming et al. 2008). The medium-mass M star 
disk does not form any giant planets at all, consistent 
with the observed paucity of gas giants orbiting M dwarfs 
(Johnson et al. 2007). 

We note that critical core mass is not constant through- 
out the disk, but instead depends on the isolation mass 
of planetesimals in the feeding zone. In the outer disk, 
where the isolation mass can reach lOOM^ or more, the 
solid accretion rate remains high throughout planet for- 
mation and the kinetic energy deposited by planetesi- 
mals inhibits the protoplanet atmosphere's contraction. 
By using lOM^ as the critical core mass throughout our 
experiments, we are overestimating the chances of giant 
planet formation in the outer disk: simulations of planet 
formation in the trans-Saturnian solar nebula show that 
critical mass can exceed 40M ffi (Dodson-Robinson et al. 
2008; Dodson-Robinson & Bodenheimer 2009). 

Now, being extremely generous not only with the feed- 



ing zone setup and disk lifetime, we investigate the 
core formation efficiency in the maximum-mass nebu- 
lae. Since core accretion requires a stable, quiescent 
disk, we cannot possibly load the disk with more plan- 
etesimal mass than contained in the Q = 1.5 nebulae 
(except through second-order effects such as high metal- 
licity) . If giant planet cores cannot reach critical mass in 
our maximum-mass nebulae, they have no chance to do 
so in any physically realistic disk. Figure [2] shows that 
the maximum radius where core accretion is effective is 
35 AU in the A star disk, 31 AU in the G star disk, and 
23 AU in the M star disk. Core accretion might be re- 
sponsible for HR 8799 d, on a ~ 23 AU orbit (Marois et 
al. 2008; Fabrycky k Murray-Clay 2008), but the other 
direct imaging discoveries HR 8799 b and c and Fomal- 
haut b are well outside the core accretion planet-forming 
zone. 

Is it possible that a super- Jupiter could form on a 
~ 35 AU orbit and migrate outward to 50 AU or greater? 
Perhaps: there are two robust outward migration mech- 
anisms, planetesimal scattering and Type III migration. 
In addition, Morbidelli & Crida (2007) propose that two 
giant planets in a 2:3 mean motion resonance that open 
overlapping gaps in the disk may halt or reverse the nor- 
mally inward Type II migration process. 

Hahn & Malhotra (1999), followed by Gomes et al. 
(2004) and Tsiganis et al. (2005), found that a 30-50M® 
planetesimal disk could expand Neptune's orbit from 
23 AU to 30 AU. However, a planet must interact with 
at least its own mass in planetesimals to change its 
semimajor axis appreciably. Even our maximum-mass 
A star nebula contains < 2Mj up of planetesimals be- 
tween 35 and 100 AU. Since HR 8799 b and c are each 
at least 7Mj up , we can rule out core accretion followed 
by planetcsimal-driven migration. 

Type III migration, which is driven by torques at the 
planet's corotation resonance and may proceed inward 
or outward, is the least understood type of large-scale 
planet motion. However, Peplinski et al. (2008) show 
that outward Type III migration (a) requires a sharp 
inner disk edge located near the planet and (b) always 
reverses direction and becomes inward migration after 
the planet's semimajor axis has approximately doubled. 
Based on their work, we do not consider Type III mi- 
gration from the inner disk a viable explanation for the 
origins of HR 8799 b and Fomalhaut b. 

Our first numerical experiment demonstrates that even 
under the most favorable conditions possible — massive, 
long-lived disk; low-velocity planetesimals; a seed core 
of 0.1M® to start the accretion — gas giants cannot form 
by in situ core accretion beyond 35 AU. For lower-mass 
G and M stars, the core accretion scenario is even more 
unfavorable. The work of Peplinski et al. (2008) indi- 
cates that outward migration from the inner disk is also 
implausible. The next two experiments investigate other 
ways of placing super-Jovian planets on extremely wide 
orbits. 

4. EXPERIMENT 2: SCATTERING 

Next we consider the possibility that planetary sys- 
tems such as HR 8799 may in fact have formed via core 
accretion — but as a system that was initially more com- 
pact than currently observed and containing at least 
one additional massive planet. (Here we use the term 
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"planet" despite the significant uncertainties in mass and 
formation mechanism). The unseen, massive planet scat- 
ters the discovered planets, forming in the inner nebula 
according to the canonical core accretion theory, out to 
their current wide orbits. Since HR 8799 is the only 
multi-planet system that includes planets on ultra-wide 
orbits, we use it as a test case for this experiment. 

Our basic approach is to perform many N-body inte- 
grations of systems with a 1.5M Q A star, to match the 
mass of HR 8799, and four planets. Our goal is to de- 
termine if planet scattering can lead to systems qualita- 
tively similar to the HR 8799 system. Three of the plan- 
ets correspond to the discovered objects, with masses of 
7Mj up , 7Mj up , and 5Mj up . These masses are intention- 
ally chosen to be near the lower limit of current mass 
estimates (Marois et al. 2008), as previous studies have 
shown that systems containing lower planet masses are 
less likely to be rapidly disrupted by strong planet-planet 
scattering (Fabrycky & Murray-Clay 2008; Gozdziewski 

6 Migaszewski 2009). The putative fourth companion 
would need to have a mass equal to or larger than the 
other planets so that it can effectively scatter the other 
planets outward (Ford & Rasio 2008). Since the total 
mass in planets is already quite large compared to other 
known planetary systems (Butler et al. 2006), we set the 
mass of the putative fourth body to be 7Mj up . We in- 
tegrate each initial configuration of planets for 200 Myr, 
using the Mercury hybrid symplectic integrator (Cham- 
bers 1999). 

The initial semi-major axis for the unseen inner planet 
is drawn from a uniform distribution between 3 and 

7 AU. The lower limit is chosen in order to place the 
innermost planet as close to the star as it could plausi- 
bly have formed via core accretion. As shown in Figure 
[2] (black curve), for A stars there is a rapid rise in the 
maximum core mass near 3 AU in both the medium- and 
maximum-mass nebulae. The upper limit for the initial 
semi-major axis of the inner planet is based on the need 
for the outer planets to be close enough to the central star 
for core accretion to remain a viable formation mecha- 
nism even when the planets are placed several mutual 
Hill radii apart (see below). 

We assign initial semi-major axes for the three outer 
planets such that the difference in the semimajor axis 
of each each pair of neighboring planets is a^+i — = 
K x TH,i,i+i, where K is a constant of order unity and 
TH,i,i+i is the mutual Hill radius of the «th and i + 1th 
planet (Chatterjee et al. 2008): 
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The choice of K affects the spacing of the planets and 
thus the timescale until the first close encounter. Based 
on similar scattering experiments (but with three less 
massive planets), we expect that a small choice for K 
(e.g., ~ 2 — 3.5) is likely to result in a very rapid first 
scattering event, which would cast doubt on the plausibil- 
ity on the initial conditions. A somewhat larger choice of 
K r*> 4 allows the system to evolve for hundreds or thou- 
sands of orbital periods before the first close encounter. 
Thus, the scattering is triggered by dynamical instability 
in a chaotic system rather than the choice of initial con- 
ditions (Chatterjee et al. 2008). In this case, we cannot 



incrase K much further, as the initial semi-major axis of 
the outermost planet becomes prohibitively large in the 
context of formation via core accretion (Figure 2). 

Each planet is assigned an initial eccentricity of less 
than 0.05, an initial inclination of less than 0.025 radi- 
ans, and random angles for the initial pericenter direc- 
tion, initial longitude of ascending node and initial mean 
anomaly. We perform two sets of ~300 simulations, one 
with K = 4.0 and another with K — 4.5. The sequence 
of the three observed planets is chosen randomly. Each 
planet is assigned an initial eccentricity of less than 0.05, 
an initial inclination of less than 0.025 radians, and ran- 
dom angles for the initial pericenter direction, initial lon- 
gitude of ascending node and initial mean anomaly. 

The median initial semi-major axis of the outer-most 
planet is 33.8 ± 5.0 AU for K = 4.0 and 40.2 ± 7.0 AU 
for K = 4.5. Thus, three of the simulated planets typi- 
cally start at smaller semi-major axes than the observed 
separations for all three observed companions. One sim- 
ulated planet typically has an initial semi-major axis be- 
tween HR 8799 c and d in the K — 4.0 simulations, or 
comparable to the present separation of planet c in the 
K — 4.5 simulations. As shown in Figure^] (black curve), 
the maximum core mass cannot reach critical core mass 
beyond separations of ~ 35 AU. If taken literally, this 
would suggest that only our simulations with an inner 
planet initially near 3 AU could result in four giant plan- 
ets with each separated from its neighbors by at least 
4 mutual Hill radii. Nevertheless, we consider initial 
conditions that include systems with somewhat wider 
initial spacings, given the large current separations of 
planets HR8799 b, c and d. While the initial conditions 
with K = 4.0 are more consistent with the modeling in 
most of these simulations result in rapid instabilities 
(within ^lMyr or ^5,000 orbital periods of the outer 
planet). This motivated the additional simulations with 
K = 4.5 which had their first close encounter on a longer 
timescale of ~24Myr, near the minimum age of the host 
star (30 Myr; Marois et al. 2008). 

For the simulations with K = 4.5, ~ 32% of the sys- 
tems retained exactly three planets at the end of the inte- 
grations. With K = 4.0, the retention rate was far worse: 
< 1% of the systems retained exactly three planets. In 
both cases, the systems that retained three planets did 
not have architectures similar to HR 8799, since one of 
the planets had migrated inwards and thus only two re- 
mained at large separations. Some of these systems did 
not suffer an orbital instability (so no planets at separa- 
tions comparable to HR 8799b), while others had planets 
with final semi-major axes quite different from those of 
HR 8799 (e.g., outermost planet with a semi-major axis 
much greater than the other planets). Thus, our simula- 
tions suggest that planet scattering alone cannot produce 
a stable system that looks like HR 8799 — although Fo- 
malhaut b, with a solitary planet on a wide orbit, might 
be possible. 

The above result leads us to consider the possibility 
that HR 8799, and future discoveries of multiple, mas- 
sive planets on wide orbits, might be dynamically un- 
stable. Given the relatively young age of the host star 
(< 160 Myr; Marois et al. 2008), we might be observing 
the planetary system in a transient phase while planet- 
planet interactions are still acting to resolve a recent 
dynamical instability (Gozdziewsi & Migaszewski 2009). 
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Thus, we analyzed snapshots of each system taken every 
200,000 years up to a maximum duration of 160 Myr. 
We include snapshots at times as small as 200,000 year 
after the beginning of our n-body integrations. 

We identify those snapshots with three outer planets 
that have semi-major axes between 20 and 100 AU, qual- 
itatively similar to HR 8799 b, c and d, plus a fourth 
planet that is closer to the star and would have eluded 
detection. For the simulations with K = 4.0 (4.5), we 
find - 0.01% - 0.07% (~ 0.01% - 0.02%) of snapshots 
meet this criteria depending on whether we look at the 
entire 160 Myr or just the first 30 Myr of our simulations. 
These simulations suggest that if the HR 8799 system is 
the result of core accretion at smaller separations fol- 
lowed by planet-planet scattering, then we would be ob- 
serving the system during a very brief transient phase 
of its evolution. Thus, planet scattering is not a viable 
formation model if systems such as HR 8799 are found 
to be common. However, at least for the time being, HR 
8799 is the only multi-planet system discovered by direct 
imaging and likely subject to strong selection effects. 

In summary, we cannot conclusively rule out the model 
of core accretion followed by planet scattering for young 
systems that might still be dynamically evolving. Our 
results are consistent with the findings of Gozdziewski & 
Migaszewski (2009), Scharf & Menou (2009) and Veras 
et al. (2009), who find dynamical relaxation timescales of 
only 10 Myr for tightly packed systems. Indeed, Scharf 
& Menou report a planet retention rate of 20% after 
10 Myr. Furthermore, Raymond et al. (2009) show that 
while planetesimals may help retain planets that would 
otherwise have been ejected, the eccentricity distribution 
of massive scattered planets remains wide — contrary to 
the case of HR 8799 where stability requires nearly cir- 
cular orbits (Lafreniere et al. 2009). Even Fomalhaut 
b, with a minimum eccentricity of 0.13 (Kalas et al. 
2008), cannot have a periapse-apoapse swing of more 
than 15 AU without disrupting its neighboring debris 
disk. We recommend that observers examine the planet 
detection rate as a function of stellar age, controlling for 
the planets' dimming with time. If planet detection rate 
is found to be independent of stellar age, it would confirm 
our finding that planet scattering is not the dominant 
mode of producing detectable planets on wide orbits. 

Finally, there is the possibility — not explored here — 
that large planet masses may allow for significant planet- 
planet interactions to begin occurring while there are still 
significant planet-disk interactions (e.g., migration, ec- 
centricity damping, and accretion) that are not included 
in our model. Unfortunately, detailed models of planet- 
disk interactions over ~100 Myr would be extremely 
computationally demanding and are thus left for a fu- 
ture study. 

5. EXPERIMENT 3: GRAVITATIONAL INSTABILITY 

Our planet scattering experiment was not successful 
at making a stable system with massive planets on wide 
orbits, and we showed in ^3] that core accretion plus mi- 
gration have severe difficulties producing massive gas gi- 
ants on wide orbits. In this section we investigate the 
one remaining possibility, gravitational instability in the 
disk. Following (1) the breakup of a disk into fragments; 
(2) the orbital evolution of each fragment; (3) the frag- 
ments' late-stage accretion of both planetesimals and gas; 



and (4) the contraction of a fragment to approximately 
Jupiter's radius is beyond the capability of any single nu- 
merical model. As a first step we simply verify that our 
model disks can become gravitationally unstable. 

5.1. Mathematical Method 

We search for exponentially growing spiral modes us- 
ing the linear stability analysis method of Adams et 
al. (1989), further developed by Laughlin & Rozyczka 
(1996). The method is described extensively in the liter- 
ature, so we provide only an overview here. 

The disk surface density is decomposed into an unper- 
turbed component plus an exponentially growing mode 
in spiral form: 

E(r, 9, t) = E (r) + Ei(r) exp[i(ujt - m6)]. (12) 

In Equation ll2l So is the unperturped surface density, Si 
is the perturbed component, m is the number of spiral 
arms, and (r, 9, t) are the radial, azimuthal and time 
coordinates. The complex eigenfrequency lo is given by 



P - H, 



(13) 



where Q p is the fixed pattern speed and 7 is the growth 
rate of the spiral wave. The disk radial velocity, az- 
imuthal velocity and potential ip due to self-gravity 
can be modally decomposed in an analogous manner to 
the surface density. Unperturbed quantities So and c 2 s 
(sound speed) are specified by the disk models described 
in <J21 Here we assume that the disk is vertically isother- 
mal, consistent with the Chiang & Goldreich (1997) pas- 
sive disk model adopted in EJ3J 

We begin our search for growing modes by calculating 
the rotation curve — including self-gravity and pressure 
support — and corresponding epicyclic frequency: 



rVl 2 {r) 



GM* 



c 2 dS d 
— H 

Sq dr dr 



(V>o) 



K 2( r ) = 1A. 



(r 2 ny 



(14) 



The disk self-gravity term, the third term on the right- 
hand side of Eq. [141 is numerically evaluated as 



1 de 



GZ (r')r'dr' 



1/2" 

'R\n (r 2 + r' 2 — 2rr'cos6 + rj 2 ) ' 

(15) 

where the softening parameter 77 = 0.1 allows the az- 
imuthal integral, which diverges at the disk edges, to be 
evaluated. 

The next step is to calculate the potential due to self- 
gravity of the growing spiral perturbation: 



■Rout fit! 



ipi = -G I Ei(r')r'dr' / d ° . 

Jr Jo \Jr 2 + r' 2 — 2rr' cos 9 

(16) 

We define the perturbed component of the disk enthalpy, 
which works against the growth of the spiral mode, as 



2 Si 



hi = c s — . 
^0 



(17) 



The linear stability analysis method is built on the ap- 
proximation that density waves propagate only radially 
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and azimuthally. Such a formulation requires that we 
treat the disk as a two-dimensional surface, so we can- 
not use volume density in our calculations. Therefore, 
although we have introduced enthalpy into the calcula- 
tion to treat adiabatic density perturbations, we use the 
isothermal sound speed, which does not depend on den- 
sity, in Equations [T4l and [T71 

Finally, we combine the linearized equations of motion 
(continuity, radial force and azimuthal force; written out 
in full in Adams et al. 1989) into a single homogeneous 
integro-differential equation in terms of hi and tp\: 

d 2 d 

—z (hi + i/ji)+A— (hi + i>i) + B(hi + ipi) + Chi = 0, 
dr z dr 

(18) 

with the coefficients A, B, C and v defined as 



A=-\o. 
dr 



B 



C=- 



m 2 Am Q, i 



dv 2m £1 d f k 
log 



r 2 r 2 k (1 — v 2 ) dr rv k dr 

K 2 {l-V 2 ) 



u) — mf2 



(19) 

For the full derivation of the governing equation, see 
Adams et al. (1989). 

The goal of linear stability analysis is to find combina- 
tions of m, ttp and 7 that satisfy Equation llSl Such com- 
binations represent exponentially growing spiral modes. 
Since we do not a priori know the perturbed surface den- 
sity Ei associated with any mode, we discretize the disk 
and write the governing equation in matrix form: 



M jk E lh = 0. 



(20) 



In Equation [201 Ei fc is the perturbed surface density at 
position A: of a discrete, linear grid. The matrix Mjk acts 
on the perturbed surface density to evaluate the left-hand 
side of Equation [TBI Since Equation [15] is homogeneous, 
it has a solution if the complex eigenfrequency lj (Equa- 
tion [T2]) is an eigenvalue of Mjk, such that 



det\M jk \ = 0. 



(21) 



To construct Mjk, we identify the rows with variable 
r in Equation llGI and the columns with variable r' . Mjk 
must encode a method of integrating Poisson's equation 
(Eq. [To]) , which is singular on the diagonal, to calculate 
•ipx. For the diagonal elements of Mjk, we use open- 
interval Romberg integration (Press et al. 1992) to eval- 
uate the azimuthal integral. We incorporate coefficients 
of the extended trapezoidal rule to perform the integra- 
tion over r' . The top and bottom rows of matrix Mjk 
encode the boundary conditions describing the flow of 
gas at the inner and outer disk edges. We follow Laugh- 
lin & Korchagin (1996) in requiring the radial component 
of the velocity to vanish at both boundaries. The effect 
of these boundary conditions is to forbid the disk from 



expanding inward or outward, which results in lossless 
reflection of density waves at the disk boundaries. The 
rigid boundary approximation and the value of the soft- 
ening parameter rj in Eq.[T5] combine to artificially affect 
stability at the disk edges (Adams et al. 1989). 

Finally, we need a way of honing in on eigenvalues that 
provide solutions to Equation [THJ We restrict our analy- 
sis to two-armed spirals and do not consider higher-order 
modes. Viable pattern speeds have their corotation res- 
onances inside the disk, so we check values of Q p that lie 
along the rotation curve (Equation[JJJ. We desire modes 
that grow within a few to a few hundred orbital peri- 
ods, so they can substantially modify the structure of the 
disk before accretion or stellar heating stabilizes it: the 
growth rate is therefore constrained by 0.01fl p < 7 < f2 p . 
We use the Newton- Raphson method to identify the valid 
eigenmodes beginning with gridded guesses for f2 p and 7. 
Once we have found a solution to Equation[3T] we set the 
first component of the perturbed surface density vector 
to unity and use least-squares analysis to solve the over- 
constrained system resulting from Equation [20l 

5.2. Results 

Our linear stability analysis reveals that the maximum- 
mass nebulae do indeed undergo global gravitationally 
instabilities. In all the maximum-mass disks, regardless 
of host star mass, we found multiple m = 2 growing 
modes with corotation resonances inside the disk. We 
list the fastest-growing modes uncovered by the Newton- 
Raphson algorithm — those with growth rates greater 
than 10% of the orbital speed at the outer disk edge — in 
Table [JJ For each mode the maximum- amplitude pertur- 
bations straddle the corotation resonance and are con- 
tained within the inner and outer Lindblad resonances, 
which occur at fl p — fi=p k/2. Figure [3] shows six of the 
modes listed in Tabled] Units on the color scale are arbi- 
trary as the modes grow exponentially. In our simple for- 
mulation, where the disk characteristics are completely 
determined by the stellar radius and temperature, the 
viability of spiral modes does not appear to depend on 
the mass of the star: any system with a high disk/star 
mass ratio can become unstable. 

Since our maximum-mass disks are both thin and mas- 
sive relative to their host stars, containing between 1/4 
(A star) and 1/3 (M star) of the total system mass, the 
fact that they become unstable is not at all surprising. 
Numerous previous investigations, including Adams et 
al. (1989), Laughlin et al. (1998), Lodato & Rice (2005), 
Kratter et al. (2008), Cai et al. 2008, and Stamatellos 
& Whitworth (2009), have verified that cold disks with 
more than 10% of the mass of their stars are subject to 
long- wavelength, global gravitational instabilities. 

Again unsurprisingly, we were not able to find growing 
modes in our medium-mass nebulae within the limita- 
tions of a grid-based search for modes using the Newton- 
Raphson scheme. If planets on wide orbits form by 
gravitational instability, they can do so only in massive 
disks. They must be relics of the early (Class 0-Class 
1) stages of protostar evolution, where massive disks are 
still present. 

At long last one of our experiments has revealed that 
there is a chance of forming the companions to HR 8799 
and Fomalhaut, and placing them on stable orbits, by a 
known mechanism. However, there is still a lot of work 
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to be done to confirm that these planets are products of 
gravitational instability. Our experiments have merely 
verified that our maximum-mass disks can become unsta- 
ble. We do not yet know what happens to the resulting 
spiral modes. We must rely on the work of previous in- 
vestigators, whose simulations have shown spiral arms 
in massive disks can fragmenting into self-gravitating 
clumps of several Jupiter masses (e.g. Boss 2000; Pickett 
et al. 2003; Mayer et al. 2004; Boley 2009). Even more 
importantly, simulations by Mayer et al. (2004) and Boss 
(2005) indicate that the protoplanetary clumps can sur- 
vive without migrating into the star. 

We caution that the viability of disk instability as a 
planet formation mechanism has not been fully estab- 
lished. The fragment survival issue requires further in- 
vestigation, as Mejfa et al. (2005) show that clumps may 
be sheared apart within an orbit or less. Pickett & 
Durisen (2007) caution that clump longevity in numerical 
simulations strongly depends on artificial viscosity, which 
is implemented differently in smoothed particle hydrody- 
namics and grid-based calculations. Boley et al. (2007) 
and Cai et al. (2009) improve upon the radiative transfer 
treatment of Boss (2007), who found rapid disk cooling 
and fragmentation, and recover a cooling rate too slow 
for any fragmentation to occur over several outer disk ro- 
tations. The emerging consensus seems to be that while 
fragmentation cannot occur in the inner disk, a < 20- 
40 AU (Rafikov 2005; Boley & Durisen 2008; Stamatellos 
& Whitworth 2008; Cai et al. 2009; Forgan et al. 2009), 
it does occur in extended disks outside of a few tens of 
AU (Boley 2009; Stamatellos & Whitworth 2009) pro- 
vided that the disks are massive. The analytical work 
of Clarke (2009) demonstrates that disks accreting from 
the protostellar core inevitably fragment upon reaching 
a radial extent of 70 AU. 

At this point we have only circumstantial evidence that 
massive planets on wide orbits, beyond 35 AU, form 
by gravitational instability. In this work disk instabil- 
ity emerges as the only possible theory from a "survival 
of the fittest" experiment, but it is possible that some 
fundamental piece of physics is missing from our un- 
derstanding of planet formation. Nevertheless, unless 
the wide-orbit planetary systems we have observed are 
dynamically unstable — unlikely given their apparently 
near-circular orbits — it seems like there is no other pos- 
sibility. Core accretion at large stellar separations un- 
equivocally does not work. Our numerical experiments 
demonstrate that gravitational instability is the only vi- 
able planet formation mechanism for planets on wide 
(a > 35 AU) orbits. We predict that HR 8799 b, c 
and d, Fomalhaut b, and similar planets that may yet 
be discovered are the result of gravitational instabilities 
in protostellar disks. 

Boley (2009) argues for two modes of planet formation, 
with core accretion dominating near the star and disk in- 
stability forming a population of gas giants at large radii. 
To Boley's conclusion that super-Jupiters on wide orbits 
should be a significant, or even dominant, component 
of the gas giant population of low-metallicity stars — in 
whose disks core accretion is inefficient — we add that the 
same is true of low-mass stars. Figure [2] demonstrates 
that gas giant formation by core accretion is all but im- 
possible for M dwarfs. Core accretion has a chance only 
in M-star disks with extremely high mass and/or metal- 



licity, but high-mass disks are likely to fragment and de- 
stroy the quiescent feeding zones necessary for steady 
growth by core accretion. 

Clarke (2009) argues that companions formed by disk 
fragmentation typically have at least 10% of the primary 
star mass. The minimum-mass disk-formed companion 
to a 0.1M© star would therefore be 10 Jupiter masses 
and would be classified as a planet rather than a brown 
dwarf. We predict that the occurrence ratio of long-period 
to short-period gas giants should be highest for M dwarfs, 
reflecting both the inefficiency of short-period planet for- 
mation by core accretion and the planetary-scale mini- 
mum mass of wide-orbit fragments formed by disk insta- 
bility. 

Having finally identified a plausible formation mecha- 
nism for massive planets in wide oribts, we conclude our 
article with a review of our findings and recommenda- 
tions for observers on how to confirm the gravitational 
instability hypothesis. 

6. CONCLUSIONS AND TESTABLE PREDICTIONS 

We have investigated the likelihood of forming mas- 
sive planets on wide orbits by three methods: near-m 
situ core accretion, scattering by an unseen companion, 
and gravitational instability. Highlighting the theoreti- 
cal difficulties presented by the companions to Fomalhaut 
and HR 8799, the prospects for forming such planets are 
generally poor. We can certainly rule out core accretion, 
as it is hard to imagine a way around the fact that long 
dynamical times destroy the chances of building critical- 
mass cores. So far, numerical simulations of Type III 
migration do not indicate that there is any possibility 
of moving a planet all the way from 35 AU — the largest 
possible stellar separation at which gas giants can form 
by core accretion — to > 60 AU (Peplihski et al. 2008). 

Planet-planet scattering can at least move a planet into 
a wide orbit; the difficulty is keeping it there. Our N- 
body simulations (Q demonstrate that once scattered, 
planets tend to leave the system. Planetesimals — not 
included in our calculations — can stabilize the orbits 
of scattered planets through dynamical friction. How- 
ever, a planet needs to interact with of order its own 
mass in planetesimals to circularize its orbit. Even our 
maximum-mass A star nebula, with total mass O.SIM© 
within 100 AU, does not contain 7-10 Mj up of planetes- 
imals beyond 35 AU. If HR 8799 is an unstable system, 
planet-planet scattering may be responsible for its ar- 
chitecture, but the shepherded debris disk surrounding 
Fomalhaut (Kalas et al. 2005) makes its present dynam- 
ical instability unlikely. 

The one planet formation mechanism that can succeed 
at creating massive gas giants on wide, near-circular or- 
bits is gravitational instability. In <j5] we demonstrated 
that the same maximum-mass disk models which could 
not form planets on wide orbits by core accretion do un- 
dergo global spiral instability. We predict that massive 
gas giants on wide orbits are the products of protostel- 
lar disk fragmentation. This prediction, however, comes 
with the caveat that the detailed physics of the disk in- 
stability process must be further analyzed and the clump 
shearing and disk cooling time problems resolved (Mejfa 
et al. [2005]; Cai et al. [2006]; Pickett & Durisen [2007]; 
Rafikov [2009]). 

A simple observational test would resolve the question 
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of how massive planets on wide orbits form. So far 11 
planets or Y dwarfs orbiting main-sequence stars have 
been discovered by direct imaging. Once the number of 
direct imaging discoveries reaches ~ 30, we can begin 
to look for simple statistical patterns in the properties 
of planet hosts. We propose dividing the stars surveyed 
by direct imaging into young and old age bins, with the 
dividing age of order ~ 100 Myr (the longest relaxation 
time in our simulations), and calculating the wide-orbit 
planet occurrence rate of young vs. old stars. If massive 
planets in wide orbits are equally common around young 
and old stars, gravitational instability is almost certainly 
their formation mechanism. Any age trend would indi- 
cate that planets on wide orbits are transient relics of 
scattering from the inner disk. In conducting such an 
experiment, care must be taken to control for the fact 
that planets cool and become less detectable over time. 

Further evidence for dual modes of planet formation, 
with core accretion operating near the star and gravita- 
tional instability dominant at large distance, would come 
from examining the occurrence ratio of long-period to 
short-period gas giants as a function of stellar mass. We 
predict that this ratio should be highest for M dwarfs, 
for which core accretion is inefficient, short-period plan- 
ets are rare (Johnson et al. 2007; Cumming 2008), and 
disk fragmentation appears to produce objects of plane- 
tary mass (Clarke 2009). Note, however, that preferred 
fragment mass in extended disks is still a subject of con- 
siderable debate (e.g. Boley et al. 2009, Kratter et al. 
2009). Core accretion on its own would produce a planet 
population where the occurrence ratio of long-period to 
short-period gas giants has the opposite trend, increasing 
with stellar mass. 

As further evidence for gravitational instability, we 
note that HR 8799 is a A Bootis star, a metal-poor Pop I 
star that probably accreted a thin layer of low-metallicity 
gas after formation (Marois et al. 2008). Even if the pro- 
toplanetary disk surrounding HR 8799 was not originally 
massive, passing through a nearby cloud could have trig- 
gered Bondi-Hoyle accretion onto the disk and destabi- 
lized it (Throop k Bally 2008). 

One further discriminant of planets formed by disk in- 
stability and core accretion/scattering is near-infrared 
color. Fortney et al. (2008) find that model planets 
with the ~ 5 x Solar metallicity of Jupiter and Saturn — 
generally interpreted as a signature of core accretion — 
are redder by 1.5 mag in H-K than solar-metallicity mod- 



els. Note, however, the simulations by Helled & Schubert 
(2008) showing that planets formed by disk instability 
may later be enriched by planetesimal capture. 

Finally, direct detection of planets is still a young 
enough field that finding as many planets as possible, 
without regard to the properties of their host stars, is a 
worthwhile goal. Gravitational instability naturally pro- 
duces planets on wide orbits (20-150 AU, Mayer et al. 
2004; Boley 2009; Clarke 2009; Rafikov 2009; Stamatel- 
los & Whitworth 2009), but the core accretion zone of 
proto-A stars extends to 35 AU and can be directly im- 
aged for nearby targets. Our work confirms the finding of 
Kennedy & Kenyon (2008) that A stars have wider planet 
formation zones than lower-mass stars. Lovis & Mayor 
(2007) also found the highest average planet mass around 
A stars vs. G and M stars, which aids planet detectabil- 
ity. Furthermore, A star disks experience the highest 
infall rates during formation and are better candidates 
than lower-mass G or M stars for triggered gravitational 
instability (Kratter et al. 2009). We do not believe it is 
a coincidence that HR 8799 and Fomalhaut are both A 
stars. Where the goal is simply to image planets, A stars 
make the best targets. Note that here we are distinguish- 
ing between the bulk planet detection likelihood — which 
so far is highest for A stars in both radial velocity and 
direct imaging surveys — and the ratio of long-period to 
short-period planet occurrence. 
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TABLE 1 

Fastest-Growing Modes in Maximum-Mass Disks 



Star Type 


H at Disk Edge (s" 1 ) 


Pattern Speed a 


Growth Rate 


Corotation (AU) 


Position in Figure 3 


A 


2.93 x ltr 10 


15.5 
2.1 
5.1 


3.2 
1.1 
0.5 


15.2 
59.4 
32.3 


not pictured 
lower left 
upper left 


G 


2.44 x lCr 10 


3.9 
3.2 
4.9 


7.1 

2.9 
2.5 


38.8 
44.4 
51.6 


upper middle 
not pictured 
lower middle 


M 


1.79 x icr 1() 


5.4 
2.4 


4.2 
2.5 


30.2 
53.7 


upper right 
lower right 



a Pattern speed and growth rate arc given in units of the orbital angular speed Q at the outer disk edge. 



12 Dodson-Robinson et al. 




Fig. 1. — A passive disk in radiative equilibrium with its parent star forms the basis of our numerical experiments. Top: midplane 
temperature as a function of distance from the star for our three model protostars, the 1.5Mq A star, the IMq G star and the 0.5Mq M 
star. Bottom: Solid surface density as a function of distance. The two sets of curves represent the maximum-mass nebulae, with Toomrc 
Q = 1.5, and the medium-mass nebulae, with Q = 8. The two discontinuities in the surface density distribution are the water ice line (< 5 
AU) and the methane condensation front. As long as the nebula is optically thick, midplane temperature is independent of solid surface 
density. 
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Maximum-Mass Nebula (Q=1 .5) 



_ 1000.0 



CO 

co 100.0 



o 
O 
■< — ' 

CD 
C 
TO 

a. 



X 
CO 




10.0 



20 



40 60 
Distance (AU) 



80 



100 



CO 
CO 
CO 



01 

i 

o 
O 

ts 

c 

0_ 

E 
E 

"§ 



100.0 



10.0 



Medium-Mass Nebula (Q=8) 
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Fig. 2. — Even our model A star hosting its maximum-mass nebula cannot form gas giant planets by core accretion at distances larger 
than 35 AU. Top: Maximum planet core mass as a function of distance in the maximum-mass nebulae. Assuming a critical core mass of 
lOMgj, the giant planet-forming region is restricted to a < 35 AU in the A star disk, a < 31 AU in the G star disk, and a < 23 AU in 
the M star disk. Bottom: In the more typical medium-mass nebulae, giant planet formation is restricted to a narrow region centered at 
~ 10 AU in the model A and G star disks. The M star disk is not capable of forming any gas giants by core accretion. 
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Fig. 3. — The maximum-mass nebulae are unstable to exponentially growing m = 2 spiral waves even though they are locally stable 
to axisymmetric perturbations. The left column shows growing eigenmodes in the A star maximum-mass nebula, the middle column 
shows growing modes in the G star nebula, and the right column shows spiral structure in the M star nebula. The color scale proceeds 
smoothly from red in the overdense regions, to green in the unperturbed regions, to blue in the underdense areas. Maximum surface 
density perturbations straddle the corotation resonances, located between 15 and 53 AU from the star for the instabilities shown here. 
Our numerical experiments demonstrate that gravitational instability is the only viable planet formation mechanism for planets on wide 
(a > 35 AU) orbits. 



